Entanglement and Entropy Engineering of Atomic Two-Qubit States 
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We propose a scheme employing quantum-reservoir engineering to controllably entangle the inter- 
nal states of two atoms trapped in a high finesse optical cavity. Using laser and cavity fields to drive 
two separate Raman transitions between metastable atomic ground states, a system is realized cor- 
responding to a pair of two-state atoms coupled collectively to a squeezed reservoir. Phase-sensitive 
reservoir correlations lead to entanglement between the atoms, and, via local unitary transforma- 
tions and adjustment of the degree and purity of squeezing, one can prepare entangled mixed states 
with any allowed combination of linear entropy and entanglement of formation. 

PACS numbers: 03.65.Ud, 03.67.-a, 42.50.-p 



The properties of entangled mixed states and schemes 
for their controlled preparation are presently under vig- 
orous investigation, primarily because of their relevance 
to understanding the role of purity and entanglement in 
quantum information protocols such as quantum compu- 
tation and quantum communication jlj]. The purity and 
degree of entanglement of two-qubit states can be quan- 
tified, respectively, by the linear entropy and either the 
entanglement of formation or the concurrence {§. Here, 
we propose a scheme using interactions in cavity quantum 
electrodynamics (cavity QED) which enables the prepa- 
ration of states of two atomic qubits with any allowed 
combination of linear entropy and concurrence. 

Our scheme uses the technique of quantum-reservoir 
engineering || in a cavity QED setting to couple a pair 
of two-state atoms collectively to an effective squeezed 
reservoir. The phase-sensitive quantum correlations of 
the reservoir are transferred to the two-atom system to 
produce entangled atomic states H, The degrees of 
purity and entanglement of the atomic states can be 
controlled through the excitation time, through proper- 
ties of the effective squeezing (i.e., the degree and pu- 
rity of squeezing), and through adjustment of the rela- 
tive strengths of amplitude and phase coupling to the 
reservoir. We are thus able to scan the entire allowed 
region of the linear entropy-concurrence plane (or lin- 
ear entropy-tangle plane 0), including the region be- 
tween the Werner states jpf and the recently charac- 
terized maximally-entangled mixed states (MEMS) 
Other cavity- QED-based schemes for entangling a pair 
of atoms have been proposed and even implemented (see, 
e.g., [0, ||, [|), but these schemes have focussed primarily 
on generating maximally-entangled pure states. 

In our proposal, two atoms are assumed to be tightly 
confined inside a high-finesse optical cavity and separated 
by a distance that is sufficiently large that they can be 
individually addressed by probe lasers and so that there 
is no direct dipole-dipole interaction between them. The 
cavity has a field decay rate of n and a frequency o>, and 
may, if desired, be driven with broadband thermal light 
characterized by a mean photon number n. Two stable 
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FIG. 1: Atomic level scheme for each atom. The excited 
states have energies hujj (j = r, s, t). 



ground states ( 1 0) , 1 1) ) of each atom constitute the qubit 
states (Fig. [I]). The cavity field and two auxiliary laser 
fields drive two separate resonant Raman transitions be- 
tween these states. In particular, transitions |1) <-* |r) 
and |0) <-> \s) are driven by detuned laser fields with 
(real) Rabi frequencies Q r and Q s and relative phase dif- 
ference ip, while the transitions |0) <-> |r) and |1) <-> \s) 
are strongly coupled to the cavity mode, with coupling 
strengths g r and g s (assumed the same for both atoms). 
Detunings of the fields from the excited states \r) and 
|s) are given by A r and A s . A fifth state \t) is virtually 
excited from |0) by another strongly detuned laser field, 
adding an additional ac-Stark shift to the state |0). 

The master equation for the total system density op- 
erator is (taking h = 1) 

P = -i [H, p] + C cav p + £ sp0 nP, (1) 
where H = H cav + H &t + H int , with # cav = Lua^a, 

H ai = ^ { u r \n) (n \ + ui s \si) (si\ + u> t \U) (U\ 

i=l,2 

+ S \U) (U\ + {(n r /2)e-^rt | r .j + jj c j 
+ [(O s /2)e- i ^. t+ ^|s i )(0 i |+H.c] 
+ [(n t /2)e- <Wi t*|t i )(0 i |+H.c] }, 
H^t = {9r\ri)(0i\a + g s \si){li\a + Ti.c.), (2) 

i=l,2 
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(H.c. denotes Hermitian conjugate) and 

£cavP — k(1 + n) (2apa) — a) ap — pa' a) 
+«n (2a) pa — aa' p — paa^) . 



(3) 



Here, u>Lj (j = f,s,t) denote the laser frequencies. The 
term C sp0 nP describes atomic spontaneous emission. 

To isolate the essential dynamics, we assume large de- 
tunings of the light fields from the excited atomic states 
(i.e., |Aj| ^> Clj, g r , g 8 , K,^j, where jj is the linewidth of 
state \j)), so that atomic spontaneous emission is negli- 
gible and the excited states can be adiabatically elimi- 
nated from the problem. This leads to a reduced master 
equation for a pair of effective two-level atoms (involv- 
ing states |0) and |1)) coupled to the cavity mode. This 
reduced system is characterized by the parameters 
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where j3 r and (3 S are the two (Raman) coupling strengths, 
and rj r and w s are the ac-Stark shifts per cavity photon 
induced in |0) and |1), respectively. 

To further reduce the model, we assume the "bad- 
cavity" limit, k ^> \fl r ,s\ , \Vr,s\- This enables us to adia- 
batically eliminate the cavity mode, which yields a mas- 
ter equation for the atomic density matrix in the form 

p = (2[3 2 /k)(N+1)(2S P S^ -S^Sp-pS^S) 

+ {2{3 2 /k)N (2S^ P S - SS^p- pSS" 1 ) 

- {2(3 2 /k)M (2S t pS' t - S*S*p- p5 f 5 f ) 

- (2(3 2 /k)M* (2S P S - SSp- pSS) 

+ {rf/2K)n{n + 1) (2PpP ] - P^Pp - pP^ P) . (5) 

Here, 1 = 1 - (3 2 s , rj 2 = (r] r - r/ s ) 2 , and 
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(n+l)p 2 + nf3 2 
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(2n + 1) P r f3 s e iv 
J 2 : 



(6) 



while S = (cf + cvT) j\[2 and P = o^a^[ + a^o^ 
collective atomic operators, with o~J — |0j) 

The derivation of (||) also requires that the phase of the 
effective two-level system remains constant with respect 
to the laser phase difference ip. That is, the effective 
atomic system and squeezed reservoir must be "resonant" 
with cachothcr, which requires that 
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(7) 



It is to satisfy this condition while retaining flexibility in 
our choices of f2 riS and A r;S that we use the additional 
transition |0) «-> \t). The level shift Of/(4At) provides 
an extra degree of freedom with which to satisfy (0). 

In (||) , the terms proportional to /3 2 describe the collec- 
tive (amplitude) coupling of our pair of effective two-level 



atoms to an effective squeezed reservoir, with the de- 
gree and purity of squeezing characterized by the param- 
eters {N, M} j|, ||. In particular, the effective squeezed 
quadrature variance is proportional to (N — \M\ + 1/2) 
and ideal squeezing corresponds to \M | 2 = N(N + 1), 
which requires that n — 0. The last line of (||) describes 
phase damping of the atomic qubits caused by coherent 
scattering of off-resonant (thermal) intracavity photons. 
There is no phase damping if n = or if rj r = rj s . 

A feature of the present system is that the strengths 
of the amplitude and phase damping terms are indepen- 
dently adjustable, so that, for example, one can be made 
to dominate the other (remembering that @) must re- 
main satisfied). Also, by switching off all sources of light 
(i.e., setting /3 r = [3 S = and h — 0) the state of the two- 
atom system can in principle be "frozen" at any instant. 

To begin our analysis of we note first that as- 
sociated with the collective coupling of the atoms to 
the reservoir are certain decoherence-free states, which 
decouple completely from the dynamics |9). In par- 
ticular, defining = (|00) ± |11» /\f2 and |i/; ± ) = 
(|01) ± 1 10)) /V2, one finds that decouples for all 
parameter choices, while decouples if N = M = 0. 

As a first example, we consider the case in which phase 
damping can be neglected (i.e., rj — or rj 2 <C (3 2 ). The 
steady state density matrix p ss is then, assuming an ini- 
tial state that has no projection onto given by 



( pii pi4 \ 

p 2 2 P23 
PSS ~ p 32 p 33 ' 
\ P41 p 44 / 

specified in the basis {|11) , 1 10) , |01) , |00)}, with 
\M\ 2 (1 - 2N) + N 2 (l + 2N) 



(8) 



Pn = 



(l + 2N)L 



1 1 



P22-P33-P23---—, PU~ {1 + m)L 



(9) 

(10) 



where L = l + 3N(l + N) - 3\M\ 2 . 

Examples of the steady state and time evolution to 
the steady state are shown in Figs. || and S, respectively, 
plotted as points in the linear entropy-concurrence plane 
pp| . We plot a minor variation of the true definition of 
the concurrence and call it the free concurrence Cf rcc = 
(Ai — A 2 — A3 — A 4 ), where Ai_ 4 are the square roots of 
the eigenvalues, in decreasing order, of pp, where p = 
(o y ® <J y )p* (cry ® a y ) , with o~ y — — i{a + — <7 _ ) . The use of 
the free concurrence enables separable states (Cf rco < 0) 
to be more readily distinguished. The linear entropy is 
given by S'l(p) = (4/3)[l — Tr(p 2 )]. On each graph we 
also plot lines corresponding to the Werner states, pw — 
£ (0+| + (1/4) (1-014 (0 < £ < 1), the MEMS of [§, 
which have the maximum amount of entanglement for a 
given linear entropy, and thermal states p t h = [C |0) (0| + 
(1 - C) |1) (1|]® 2 (0 < C < !)■ For ideal squeezing the 
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FIG. 2: Steady state values of Cf rec (p) and Sl(p) for selected 
values of N and < \M\ 2 < N(N + 1), for initial state |00). 
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FIG. 3: Evolution of Ct ree (p) and Sl(p) to the steady state 
with ideal squeezing (|M| 2 = N(N + 1)) for the following 
initial states and values of N: { 1 00) ,N = 0.2} (o), {|00} , N = 
0.5} (□), {\00},N = 1} (o), {|00),7V = 5} (•), {|11),JV = 
0.2} (A), {|01) , TV = 2} (x), {|01> , AT = 0.01} (+). Note 
that the points on each curve are not equally spaced in time. 



steady state described by fel) is the pure state 
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(11) 



These states lie on the left-hand vertical axes of the fig- 
ures and approach the Bell states |</>+) in the limit of 
large squeezing (i.e., large N, and tp — ir or 0). Nonideal 
squeezing (Fig. ||) generates steady states that can lie es- 
sentially anywhere below the Werner line. Note that for 
large N the steady states closely approximate mixtures 
of (</? = tt) and pf = diag {1/3, 1/6, 1/6, 1/3}. 

Time evolution with ideal squeezing from initial states 
with zero projection onto \ip~) can also sweep out the 
region beneath the Werner line (Fig. |J). When the ini- 
tial state does have a projection onto \ip~) (e.g., 1 1) ) , an 
interesting range of points on the plane can also be ac- 
cessed, including an area above the Werner line and the 
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FIG. 4: Evolution of Cf rcc (p) and Sl(p) for initial states ( |l^ ) 
with 9 = {tt/4, 7T-/2}, N = 3.1, and n = 0. The dotted curves 
show evolution produced by phase decay turned on after ap- 
plication of the unitary transformation U to a selection of 
states from the = tt/2 curve. For this evolution, amplitude 
coupling is disabled, but rj ^ and n = 1. 



region along the boundary at Cf reo = between separable 
and entangled (including the maximally mixed entangled 
state at the intersection of the Werner and MEMS lines). 

States above the Werner line can also be generated by 
initially preparing the separable pure superposition state 



|tf (0)> = {cos(r5/2)l 2 



ism{0/2)<7 y }® 2 \OO}. 



(12) 



with < 8 < tt/2, and then applying the effective reser- 
voir interaction with strong, ideal squeezing. In this 
case, the time evolution follows paths as shown in Fig. |^. 
There is a small region of the plane above the Werner 
line which cannot be reached by this method. However, 
this region can be accessed by switching off the squeezed 
reservoir interaction ((3 r ,s = 0) and employing phase de- 
cay (77, n 0) from initial states prepared on the 9 = tt/2 
curve of Fig. ^ and to which the local unitary transforma- 
tion U = (1/2) {a x + <j z } ® {I2 — i(Jy} (requiring single- 
atom addressing with appropriate laser Raman pulses) 
has first been applied. This is also illustrated in Fig. ||. 

We turn now to practical issues associated with our 
scheme. Analysis of (^) shows that the slowest rate fea- 
turing in the dynamics is (4/? 2 / 'k)(2N - 2\M\ + 1), which 
exhibits the characteristic inhibited decay associated 
with atomic damping by a squeezed reservoir ]Tl[ . Mean- 
while, inclusion of atomic spontaneous emission effects 
(due to finite excited state populations) into the model 
reveals characteristic rates /2A|) (J — r, s, t). Tak- 
ing the rate for j = r to be the maximum, and set- 
ting n — 0, the condition that spontaneous emission be 
negligible during the state preparation period reduces 
to 2g2/( 7rK ) > [1 - ^/N/{N + I)]" 2 . This amounts 
to the condition of strong coupling cavity QED, made 
somewhat more stringent however owing to the inhib- 
ited atomic decay rate. If we consider a recent cav- 
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FIG. 5: Evolution of Cf rcc (p) and Sh(p) with spon- 
taneous emission effects for (g r , g s , t, jj,Cl s , Aj) /2tt = 
(110, 110, 14.2, 5.2, 100, 8000) MHz (j = r, s, t), Of = Q? r - fij, 
and n = 0, with the following initial states and values of Q r : 
{Eq.(O), 6» = tt/2, fi r = 120 (AT = 2.3)} (•), {| 00), fir- = 
173 [N = 0.5)} (o), {|ll),fi r = 245 (N = 0.2)} (A), 
{|01},fi r = 110 (N = 5)} (x), {|01),n r = 458 (N = 0.05)} 
(+). Note that the points are not equally space in time. 



concurrence can also be formulated using actual squeezed 
light to drive the cavity and just a single cavity-mediated 
Raman transition. Further, with nondegenerate squeezed 
light fields driving distinct cavities it is possible to engi- 
neer the states of distantly separated atoms Q . 

The scheme presented here can also be used to prepare 
spin-squeezed states of a larger number of cavity-confined 
atoms (^]. While preparing this manuscript, we became 
aware of proposals for preparing such states which use 
essentially the same configuration as in this work, only 
not in the guise of quantum reservoir engineering []l9f . 

In conclusion, we have proposed a scheme for engineer- 
ing atomic two-qubit states with any allowed combina- 
tion of linear entropy and concurrence. This should open 
the door to detailed experimental investigation of purity 
and entanglement in quantum information protocols. 

We acknowledge helpful discussions with S. Rebic and 
W. J. Munro. This work was supported in part by the 
Marsden Fund of the Royal Society of New Zealand. 



ity QED experiment for which the parameters achieved 
were (#,«;, 7) /2tt = (110, 14.2, 5.2) MHz @, then for 
N = 2 the above inequality reads as 332 ^> 30, indi- 
cating that sufficiently strong coupling is experimentally 
realistic for achieving significant levels of effective squeez- 
ing. Furthermore, setting, e.g., O r /A r = 0.02 and using 
the above parameters, the characteristic state prepara- 
tion time is < 50 [is, which is orders of magnitude less 
than single-atom trapping times in tightly-confining op- 
tical dipole traps (see, e.g., [jl3|, [l4|, [ll| 

In Fig. I we present sample evolutions from a model in 
which the excited atomic states have been adiabatically 
eliminated, but in which effects of spontaneous emission 
plus the cavity mode dynamics are included. Using the 
cavity QED parameters quoted above, we see that a large 
area of the linear entropy-concurrence plane can be ac- 
cessed. With the inclusion of spontaneous emission, de- 
cay into the (weakly-coupled) state can occur for 
states with no initial projection onto \ip ~). This limits 
the maximal attainable concurrence and in Fig. ^| leads 
also, for the cases with | p(0) = 0, to a very slow 
decay of Cf ree (p) and increase in S^(p) after rapid initial 
evolution to the optimal value of Cf ree (p). Note again 
though that the evolution can be frozen at any point by 
simply turning off all of the light fields. Note also that 
for the highest value of Ci ICC (p) attained in Fig. ||, one 
finds (4> + \ p \4> + ) — 0.9. Higher cavity finesses (i.e., lower 
k v alues) improve the performance of the scheme further 
and 
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